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Model calculations are presented which predict whether or not an arbi- 
trary gas experiences significant absorption within carbon nanotubes and/or 
bundles of nanotubes. The potentials used in these calculations assume a 
conventional form, based on a sum of two-body interactions with individual 
carbon atoms; the latter employ energy and distance parameters which are 
derived from empirical combining rules. The results confirm intuitive expecta- 
tion that small atoms and molecules are absorbed within both the interstitial 
channels and the tubes, while large atoms and molecules are absorbed almost 
exclusively within the tubes. 

I. INTRODUCTION 

The absorption of gases in nanopores is a subject of growing experimental and theoretical 
interest, stimulated by both fundamental scientific questions and the potential for many 
technologies [p^pH]. One of the most important questions to be addressed is whether or 
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not a specific gas is significantly absorbed within carbon nanotubes; we will define the 
word "significant" in Eq. ^ below. While the answer depends in detail on the specific 
thermodynamic conditions of the coexisting vapor (pressure P and temperature T), one 
expects that intuitive considerations based on size and energy scales ought to provide useful 
qualitative insights. For example, it has been demonstrated that gases whose condensed 



phases possess low surface tensions are strongly imbibed in these tubes |]2^ . This important 
result can be understood from either the Kelvin equation or a comparison of competing 
interaction (adhesive vs. cohesive) energies. These considerations arise in the analogous 
problem of wetting transitions ||23|| . 

This paper addresses this basic question by employing a simple, but plausible, model of 
the interaction potential, from which we compute the adsorption as a function of P and T. 
We assume that the adsorption potential can be derived from a sum of Lennard- Jones (LJ) 
two-body interactions between the host C atoms and the adsorbate. This pair potential has 
distance and energy parameters obtained with semiempirical combining rules from the LJ e 
and a parameters of the C atoms and the adsorbate p^-p7|: 

_ CTgg + CTcc f-.s 

"gc — 2 

<^gc = \/£gg^CC 

where "g" and "C" refer to the gas and C atoms, respectively. Estimates of the gas param- 
eters are given for some relevant systems in Table | |2^J2^;^^, while for C atoms we use 
o"cc = 3.4 A and e^c = 28 K pS]. These values are typical, but uncertain within a 15% 



range |^| 



This paper's approach is the following. We first choose a particular (somewhat arbitrary) 
criterion for calling the uptake "significant" . For example. Fig. |l|(a) and most of our work 
employ the criterion 

p* = page = 0.1 (2) 

where p = N/L is the one-dimensional (Id) density, with the number of adsorbed atoms 
and L the length of the tube, and p* is the corresponding dimensionless density. For gases 
of interest here, this criterion corresponds to a mean Id spacing of order 30 A. This is a 
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very low density. Although we do consider more stringent criteria elsewhere in this paper, 
the results do not differ qualitatively. The reason for the lack of sensitivity to the threshold 
is that once adsorption commences, it rises rapidly as a function of P (until the crowding 
effect of repulsive forces slows the variation of coverage with P). 

In the nanotube bundle geometry, adsorption can take place inside the tubes, in the 
interstitial channels, and on the outer surface of the nanotube bundle (Fig. Typical 
length scales for the triangular lattice of nanotubes in the bundle are: lattice constant 
17 A, nanotube radius 6.9 A, bundle diameter between 50 and 100 A, and bundle length 



~ 10 — 100 fim We will see that size is a critical variable determining uptake. Some key 
findings of this paper appear in Fig. |I|(a), which shows the uptake at a very small ratio of 
P to saturated vapor pressure (Po)- Small atoms and molecules (which typically have small 
values of egg) are strongly adsorbed within both the nanotubes and the narrow interstitial 
channels (IC's) between nanotubes. Larger particles, in contrast, do not "fit" within the 
IC's but do imbibe within the tubes. Perhaps a surprising feature of Fig. |I](a) is that a 
hypothetical gas with a very large value of egg adsorbs in neither place. This occurs because 
the relative tendency (compared with bulk condensation) of a gas to be absorbed within the 
tube at a given undersatumtion depends on the ratio of adhesive to cohesive energies. The 
geometric mean combining rule for egc implies that this ratio varies as the inverse square root 
of egg, so a large egg implies small uptake. This finding is qualitatively consistent with the 
empirical correlation between uptake and surface tension mentioned above. It also correlates 
with the physics determining wetting behavior of liquids for which the analogous comparison 
involves the same kind of interaction ratio ||23|| . 

This paper makes a number of simplifications in order to draw such general conclusions. 
Arguably the most drastic assumptions are that the nanotubes are infinite and perfect and 
that the nanotube bundles involve a unique species of tubes in a regular array (geometry 
unaffected by the adsorption). 

The outline of this paper is the following. Section || describes our model of the interac- 
tions. Section |T| presents the statistical mechanical model used in the calculations. Section 
^ reports our results. Section summarizes these and discusses open questions. 
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II. ADSORPTION POTENTIAL 



A basic assumption in our model is that the potential energy experienced by a molecule 
at position r can be evaluated by a summation of two-body interactions U (x) between the 
molecule and the carbon atoms comprising the tube: 

V{r) = Y,U{r-RO (3) 

i 

This assumption is made in the overwhelming majority of calculations of gas interactions 
with either graphite and carbon nanotubes. In the graphite case, many-body effects have 
been found to be ~ 15% corrections to ab initio pair potential sums ||3^. Hence, the 
empirical pair potential should be regarded as an effective pair potential. One might expect 
somewhat smaller many body contributions in the nanotube case because the molecule is 
somewhat farther from the nearest carbon atom and because the effective coordination 
number is larger in the nanotube case than on graphite. In contrast, the argument in 
the IC case leads to the prediction of a larger many body effect than on graphite. These 
expectations, however, might not be correct because the many-body expansion involves 



geometry-dependent competing terms of opposite signs |3J] and because the two body energy 
for the IC is typically of much larger magnitude than on a flat surface. 

Another key assumption made here is that the pair potential is isotropic and of LJ 
form: U{x) = 4e[((j/x)^^ — (cr/x)^]. There is ab initio and empirical evidence to the effect 
that anisotropy of the pair potential plays a role in adsorption potentials on graphite p5 



Nevertheless, most studies of adsorption on that surface neglect such an effect and use 
a LJ pair potential similar to what we use here. The final assumption is the use of an 
azimuthally and longitudinally averaged potential. The potential at distance r from the 
axis of the cylinder is then [0: 

V{r;R) =37reea' [-(-) /n(a;)Mn(x) - (-) Mx)M,{x)\ (4) 

where 9 = 0.38 is the surface density of C atoms and R is the radius of the cylinder. 
Here, x = r^/r^, and are the smaller (greater) of r and R. The function is 

defined as 1 for r < and {R/rY for r > R, with n a positive integer. Here we use the 
integrals 



We emphasize that each approximation introduces an error, but the quahtative trends ought 
to be rehable. At this time, the lack of high quahty ab initio calculations would seem to 
warrant this kind of approach. 

The IC potential is obtained by summing up the contribution from three nanotubes and 
azimuthally averaging the result. Figures show contour plots in the agg — egg plane 
of the reduced minimum of the adsorption potential {V^^^ = Vmm/^gg) for all of these 
sites. Inside both the tubes and in IC's there is a threshold agg value above which the 
potential becomes repulsive, corresponding to gases which are too big to fit in these restricted 
geometries; these thresholds are agg ~ 11.4 A (tubes) and 3.4 A (IC's) for nanotubes of 
radius 6.9 A studied here. Outside of the bundle, there are no such size constraints for the 
adsorbed atoms/molecules, as the adsorbate can always find a region in which the potential 
is attractive; at a fixed value of egg large systems yield larger due to their larger 

coordination number of C atoms. In all three cases the most negative values of V^^,^ occur 
for small values of egg. In the tube and external surface cases, but not the IC case, the most 
negative values of V^^^ occur for large agg (> 9 A). 



III. STATISTICAL MECHANICS 



Our interest is whether atoms are likely to go inside the tubes, in the interstitial channels, 
and on the outer surface of the nanotube bundle. This behavior is determined by the 
thermodynamic conditions (P, T) and microscopic parameters (especially agg relative to R). 
A key factor implicit here is the cohesive energy of the bulk phase of the adsorbate which 
determines a relevant pressure, i. e. saturated vapor pressure Pq. We construct a simple 
model for the low coverage regime of atoms inside nanotubes, neglecting the interactions 
between adsorbate atoms, while for atoms moving in the very confining IC's any density can 
be considered because of the mathematical simplicity resulting from the Id character of the 
system. We have discussed elsewhere the extreme quantum behavior of He at low T |n,S7|. 
In the present case we assume that classical statistical mechanics applies p8 |. 
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We now compute the chemical potential fi of the adsorbate. All of our calculations take 
the coexisting three-dimensional vapor to be an ideal gas, so that the chemical potential 
can be expressed in terms of pressure as /i = (3~^ln{[3P\^). Here (3~^ = ksT, and A = 
(27r^^/5/m)^/^ is the de Broglie thermal wavelength for particles of mass m. It is convenient 
to measure the chemical potential with respect to its value at saturation, /xq, 

A/i = /i-/io = /?-i/n(P/Po) (6) 

An analytical expression for Pq is available from computer simulation data of the Lennard- 
Jones system's liquid-vapor coexistence H, InP;; = 1.2629T* - 4.9095/T* - 0.15115/T*^ 
where Pq = PoCgg/egg and T* = ksT/egg are reduced quantities. 

Consider first the adsorption inside a single nanotube. The chemical potential of the 
ideal gas in an external potential can be expressed as a function of the number of adsorbed 
atoms, A^, and temperature: 

= (7) 

l^dr exp{~,3V(r)) ^ 

where the integral is performed over the volume of the nanotube. This is an application 
of Henry's law. Then, the chemical potential relative to its value at saturation assumes a 
simple form due to the cylindrical symmetry of the adsorption potential: 

^ ^''^"^ {2npPoJ^^dr''rexp{-pV{r))) 

Atoms in the narrow IC's are strongly confined to the vicinity of the axis so that a 
Id model is applicable and solvable for all densities. As previously discussed in the case 
of very small nanotubes ||T^, the transverse motion may be treated independently of the 



longitudinal motion and the chemical potential in this case has the form: 

fi = fij_ + /iid (9) 

where /i^ is the transverse contribution and fi^^ is the chemical potential of a Id gas. In 
general, = ea;p(— /5ej)), where {ej}j=o,i,--- transverse spectrum of individual 

atoms/molecules. At low T (/3(ei — eo) << 1), the ground state dominates the sum and 
/ix — cq. The ground state energy can be determined very accurately using the WKB method 
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40| , since the adsorption potential is well-approximated by a parabola in the vicinity of the 



IC axis 0]. Results for the potential well depths of various gases are shown in Table |. 
The Id chemical potential is obtained by integrating the Id Gibbs-Duhem relation 

where P^^ is the Id pressure. The particle density in the case of only nearest-neighbor 
interactions is given by the equation of state ||^,|i2[] 



^ ^ /o°" dz exp{-P[u{z) + zP,^]) 
dz z exp{-f3[u{z) + zP^a]) 

Here u(z) is the LJ potential describing the interactions between adsorbed atoms. The 
integration of eq. ^ leads to 



(3,., = In ('^^Wo°°^^ eM-m^) + ^Pm,o])- ^^^^ 



dz exp{-P[u{z) + zP^a\) 
Pidfi is an initial low pressure chosen such that the ideal gas limit is reproduced. The density 
dependence of the Id chemical potential is finally obtained by eliminating the Id pressure 
between eqs. |ll| and |12. 



As shown in Fig. ^, the external surface of the nanotube bundle also provid es an attrac- 
tive domain of adsorption. We have studied adsorption in the very attractive groove-like 
channel which runs parallel to the nanotube 3iXGS, clS shown in Figure 13 . Then, a pro- 
cedure similar to that employed in the case of IC's is applicable to computing the coverage. 
The contribution of the longitudinal motion to the chemical potential is determined in the 
same fashion as in the case of IC's. The ground state energy of the transverse motion can 
be estimated through a parabolic approximation for the adsorption potential at this site. 
Values of the ground state energy (-Eq^*) obtained in this fashion for the systems studied, as 
well as the well-depth of the adsorption potential (V^^^), are listed in Table |. 



IV. RESULTS 

The lines obtained by setting the coverage equal to the threshold criterion can be seen in 
Figures |l|(a)-(d). Figure |l|(a) shows this behavior in the case of p* = 0.1, for A/x* = Ap/egg = 
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—10 and T* = 1. As expected, small atoms or molecules (He, Ne, H2) fit easily inside both 
the tubes and IC's, while large molecules do not fit in the narrow IC's. Hypothetical (but 
nonexistent) atoms with agg < 2.5 A are adsorbed in the IC's only if their self-interaction 
energy (e^p) does not exceed a threshold value. The upper limit to the molecular size for 
adsorption inside the tubes can be seen in Fig. ^. Indeed, the experimental observation 
of Cgo molecules encapsulated in nanotubes is consistent with this expectation (as the 
point near agg ~ 9 A indicates). 

Including the effect of interactions does not affect our results significantly. In the frame- 



work of the gas-surface virial expansion [^ , 

N c^e^''Q^[l + pri{T)] (13) 
where Qi is the single particle canonical partition function, and 

The net effect of the virial correction is at most a 0.1% change of A/i; such a small magnitude 
is consistent with the expected behavior in the low pressure regime of interest here. 

The evolution of the diagram as a function of the adsorption criterion can be seen in Fig. 
|l](b). As the threshold density decreases (p* = 0.05 here) more systems satisfy the uptake 
criterion. Fig. |I](c) shows a similar effect on the diagram of an increase in chemical potential, 
to A/i* = —8. In both geometries, the altered criterion corresponds to more systems being 
allowed in the respective cavities. A different effect on the diagram occurs if the size of 
the nanotubes is changed, as shown in Fig. |l](d), under the thermodynamic conditions of 
Fig. |l](a). In the case of nanotubes with radius 8 A more atoms enter interstitial channels 
because of the larger channel space, while fewer atoms go inside the tubes because the 
adhesive energy decreases. The trends seen in Fig. |l](a)-(d) are qualitatively consistent with 
the behavior of presented in Figures |] and ^. 

In the Henry's law regime of low coverage there is a convenient way to characterize the 
variation of uptake with geometry. We compute the ratio of particle occupations in the 
nanotubes and IC's at the same P and T: 

p. X _ ^iVT Lr exp{-f3V) 

^ [egg, agg) - j^^^rexp(-/3\/) ^'^^ 
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where fNTiw) are the number of nanotubes (IC's) in the bundle and the integrations are 
over one region (assumed infinitely long). For an infinite array of nanotubes, Vmt/^ic = 1/2. 
The finiteness of the bundle changes the ratio; however, there is no qualitative effect on our 
conclusions unless the bundle is very small. This ratio depends on the two gas parameters, 
egg and agg. In order to simplify the presentation, we fit the general trend of systems in 
Table | to an empirical equation: 

^aagg + h (16) 

with values a = 147 K/A and b = 376 K. We then consider a function of one variable 

ncrj^r{elf,agg) (17) 

This ratio function is presented in Figures |^ and |^. In Fig. ^ we consider a common value 
of T* = 1, while in Fig.|^, we consider a fixed T = 77 K. The data in Fig. shows, as 
expected, that large (small) molecules adsorb preferentially in the nanotubes (IC's). Fig. |^ 
differs for small agg because at 77 K the entropic advantage of the tubes is manifested as a 
larger uptake there than is seen in Fig. ^ at the much lower T given by Eq. ([16|). 



V. CONCLUSIONS 



Model calculations were used to investigate adsorption in nanotube bundles. Simplifying 
assumptions were made, such as the pairwise summation of gas-surface interactions, the 
use of combining rules to determine energy and size parameters, and the continuum, rigid 
model of the carbon atoms of the tube. We studied mainly the regime of low coverage, 
where interactions between adsorbed atoms are omitted; in the IC case, this assumption 
was not needed, as a quasi-one dimensional approximation permits exact treatment of LJ 
interactions at finite coverages. The conclusions drawn are expected to be qualitatively 
accurate in general situations, and so they provide useful insight for experiments. The key 
result appears in Fig. indicating which molecules go where under "typical" experimental 
conditions. More general behavior can be estimated from the reduced potential curves 
(Figures ||-|D we have presented. 
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Williams and Eklund [31| have computed the H2 adsorption on the bounding surface 
of bundles containing a finite number of tubes. In some cases, this contribution can be a 
significant fraction of the total adsorption. Adsorption isotherms of classical gases on the 
external surface of the bundle is the subject of our current investigations to be reported in 
the future. 

We discuss the relevant experiments very briefly. Teizer et al. [Q studied He uptake 
and found consistency with our calculations for one-dimensional motion and the computed 
binding energy within the interstitial channels. Kuznetsova et al. |0 studied uptake of Xe 
and their data are consistent with our calculations of the uptake within the nanotubes. 
Interestingly, a recent experimental study of adsorption of methane in nanotube bundles 
concluded that significant IC adsorption occurs. This conclusion was reached from the 
fact that nanotubes were capped and the measured binding energy of CH4 determined (2570 
K) was 76% larger than that on graphite (1460 K) [^, which compared favorably with pre- 
vious estimates of the IC binding energy of II2, He and Ne [p!^-[T^. Our present calculations 
indicate, however, that the large size of CH4 prevent it from populating the narrow IC's 
significantly. In contrast, the external surface of the nanotube bundle is accessible and the 
binding energy in this case (Table |) is ~ 20% larger than the one for graphite. A more 
realistic potential exhibits corrugation, which we have neglected here; its effect is to increase 
the binding energy |[T^, but we have not undertaken that calculation as yet so no definite 
comparison is possible. 
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TABLE I. The values of the LJ parameters, egg and agg, for the gas-gas interactions, the cor- 
responding minimum (Vmin) of the adsorption potential and ground state energies (Eq) inside a 
nanotube (NT), in the interstitial channel (IC), on the external surface of the bundle (ext), and 
on a single graphite sheet (GR) is given. LJ parameters were taken from Ref. |28|, except for CH4, 
CF4 and SFe pi and Ceo P 



Gas 


^99 (K) 


^99 (A) 


'mm \-^^J 


(K) 


(K) 

'mm \-^^J 


'mm \-^^/ 


(K) 


V^^ (K) 

'mm \^^/ 


He 


10.2 


2.56 


-546 


-386 


-297 


-367 


-270 


-218 


Ne 


35.6 


2.75 


-1018 


-902 


-600 


-725 


-666 


-431 


H2 


37.0 


3.05 


-828 


-292 


-690 


-808 


-618 


-482 


Ar 


120 


3.40 


6 


228 


-1426 


-1607 


-1550 


-965 


CH4 


148 


3.45 


401 


789 


-1614 


-1809 


-1714 


-1088 


Kr 


171 


3.60 


2048 


2250 


-1836 


-2025 


-1981 


-1220 


Xe 


221 


4.10 


14786 


15054 


-2523 


-2617 


-2580 


-1593 


CF4 


157 


4.58 


36411 


36854 


-2539 


-2475 


-2433 


-1520 


SFe 


208 


5.25 


136492 


137196 


-3726 


-3307 


-3272 


-2056 


Ceo 


2300 


9.2 


52858932 


52863770 


-49071 


-21952 


-21924 


-14505 
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Figure Captions 

FIG. 1. (a) Diagram indicating regions of significant uptake (at thermodynamic conditions 
specified by Afi* = —10, T* = 1 and p* = 0.1) as a function of the adsorbate Lennard-Jones 
parameters. Gases lying in the domain deoted "TUBE" are absorbed within the nanotubes. Those 
denoted "IC" are absorbed within the interstitial channels, while those denoted "BOTH" ("NEI- 
THER") go to both places (neither place). Systems of particular interest are identified by dots 
(•), with parameters listed in Table |. 

(b) Diagram analogous to (a), except that curves shown utilize an alternative p* value for the 
threshold condition, i.e. p* = 0.05. 

(c) Same as in (a) for a different value of the chemical potential: Ap* = —8. 

(d) Same as in (a) in the case of a nanotube array with tubes of diameter 16 A. 

FIG. 2. Schematic picture of adsorption sites within and outside a nanotube bundle. For the 
external surface, the most attractive site, located at equal distance from two nanotubes, is shown 
here. Adsorbed atoms or molecules are represented by dots. 

FIG. 3. Contour plot of the reduced well-depth of the adsorption potential inside a carbon 
nanotube. The attractive isopotential curves ( — ) correspond to Vj^^^ increments of 10 starting 

from -90, while the repulsive curves ( ), from left to right, correspond to V^^^ = 20, 40, and 

80. 

FIG. 4. Same as in Fig. ^ for the interstitial channel. Repulsive curves, from left to right, 
correspond to = 10, 20, 30 and 40. 

FIG. 5. Same as in Fig. for the external surface of the nanotube bundle. 

FIG. 6. Expanded version of Fig. ||(a) showing the gas systems which absorb within a nanotube 
at p* = 0.1, Ap* = -10 and T* = 1. 
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FIG. 7. Ratio of the amount adsorbed inside a nanotube to that in an intersititial channel in 
the Henry's law (low coverage) regime. Here, T* = 1. 

FIG. 8. Same as in Fig. |, but for T = 77 K. 
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FIG. 13(a) 
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FIG. g(b) 
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FIG. 1(c) 



20 




FIG. g(d) 
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